
STOP 



Early Journal Content on JSTOR, Free to Anyone in the World 

This article is one of nearly 500,000 scholarly works digitized and made freely available to everyone in 
the world by JSTOR. 

Known as the Early Journal Content, this set of works include research articles, news, letters, and other 
writings published in more than 200 of the oldest leading academic journals. The works date from the 
mid-seventeenth to the early twentieth centuries. 

We encourage people to read and share the Early Journal Content openly and to tell others that this 
resource exists. People may post this content online or redistribute in any way for non-commercial 
purposes. 

Read more about Early Journal Content at http://about.jstor.org/participate-jstor/individuals/early- 
journal-content . 



JSTOR is a digital library of academic journals, books, and primary source objects. JSTOR helps people 
discover, use, and build upon a wide range of content through a powerful research and teaching 
platform, and preserves this content for future generations. JSTOR is part of ITHAKA, a not-for-profit 
organization that also includes Ithaka S+R and Portico. For more information about JSTOR, please 
contact support@jstor.org. 



DISCUSSION OF A METHOD FOR FINDING NUMERICAL 

SQUARE ROOTS 

By C. L. Bouton 

In a recent article * Professor E. R. Hedrick gave a method for approxi- 
mating to the positive square root of any positive number. Inasmuch as 
he merely states the method without giving proofs, and as the proof itself is a 
rather pretty exercise in the elementary theory of limits, it seems worth while 
to give a discussion of the method. The facts that the method is really 
Newton's method of approximating to the square root, and that it may easily 
be extended to the p ttl root, are of interest. Professor Hedrick's plea that 
this method be substituted for the method of extracting square root ordinarily 
used seems to be scarcely justified on the ground of rapidity of convergence, 
although there may be pedagogical advantages in its use with young pupils. 
The method has the further disadvantage that when applied to a number which 
is a perfect square, it yields the square root only as the result of an infinite 
process, unless indeed the first trial root should turn out to be the exact root ; 
whereas, the ordinary process applied to this case always terminates in a finite 
number of steps, provided that the given number is an integer, or a finite 
decimal fraction. 

The process is the following : Let a be the given number, select any posi- 
tive number n, and form in succession 

(i) / "»=i(" 1 + H- I )' 



Then we shall prove that the numbers of the sequence n l5 n„, ■ • ■, n if 
form a decreasing set, and approach ^a as their limit. The process converges 

* Western Journal of Education, vol. I, October, 1908, p. 252. 

(167) 
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rapidly enough to be useful in computing numerical square roots. For ex- 
ample, let a = 2, and choose n = 1. Then we compute in succession 

Wl = 1.5, 712 = 1.417, % = 1.414216 
Thus in three steps the error is about 2 in 10 6 . 

To prove that ' ini ^ (n^) = V«, and to estimate the error at any step, let 
us proceed as follows : First we have the identity 

whence 

(3) n { gy/o (t=l, 2, ..-)• 

That is, « lf n a , • • • are all greater than \fa, unless n was chosen equal to \Ja, 
in which case n x = n 2 = • • • = \Ja. From this point on we shall assume that 
n was not chosen equal to Va> so that in (2) and (3) the equality sign cannot 
hold. Moreover, from (1) and (3) we have 

(4) n < _ 1 -n < = wf - 1 ~ a >0 (i=2,3,...)> 

in i—\ 

so that 

(5) ««<%_! (» = 2, 3, •••)• 

Therefore the variable n t always decreases when i increases, but n 4 is al- 
ways greater than Va. Hence by a fundamental principle in the theory of 
limits * the limit of rtj exists when i = co ; let us call it Jf. We readily see 
that this limit is V«> for from (2) and (4) we have 

(6) n i -^ a = —7= (— k~ ) = TTr — (t=l, 2, ...). 

Passing to the limit for i = oo , we get 

i\r - Va = i— jL- = 0, 

JST+ \fa 

JST= Va. 
* See, for example, Osgood, Introduction to Infinite Series, p. 4. 
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This proves the convergence of the process. It is easy to set limits to 
the error at any step, and to form some idea of the rapidity of the convergence. 
We have from (6) and (3) 

(7 > 2^ < W «-V^< i ^F < Ya 

That is, if we set n i _ 1 — % = S f and «< — ^a = e t , then 

(8) £"•'* ('-!.».• ••>• 

Furthermore, from (7) and (3) 

(9) n, - yg (w ._ 1 _ n<) « l ^-rjg N', 

Let ' ._ — = r it so that r< is the relative or percentage error at the t th step. 
Then 

(10) r.Kiri^ (i=2,3, ...)• 
Repeated application of (10) shows that 

<»> ^i^p^r c- 1 . ».••■)■ 

Formula (8) shows that the actual error of any approximation is nearly 
proportional to the square of the difference between that approximation and 
the preceding one. Formula (10) gives the relation between the relative 
error at any step, and that of the preceding step. It shows that when working 
decimally the number of significant figures retained at any step should be 
double the number which are correct at the preceding step. These formulae 
both show that while the process may not converge very rapidly at first, it 
later converges with extreme rapidity. Thus in the computation for sji we 
found r s < 2 X 10~ 8 , hence r 4 < 2 X 10" 12 and r 5 < 2 X 10 ~ M . 

Formula (11) enables us to estimate the number of steps necessary to 
obtain a given degree of approximation, and furthermore shows the effect of 
the choice of the initial n. Of course the nearer n is chosen to v'« the more 
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rapid is the convergence. The method is therefore particularly useful for 
computing square roots when an approximate value is known, as for example, 
in extending a given table of square roots to a larger number of places. Thus, 
ifwetakea=3, and choose n= 1.732, a single step gives n x = 1.732050808314. 
If to this we apply the correction computed from (8) we find 1.732050807569, 
while the value of \/3 correct to thirteen decimal places is 1.7320508075689. 

Furthermore, it should be noticed that the commission of an arithmetical 
error at any stage of the work does not vitiate the ultimate result, for it 
simply amounts to starting over with a new value of n. In fact, such a mis- 
take may even increase the rapidity of the convergence. 

A geometrical proof of the convergency of the process is easily given. 
Draw the rectangular hyperbola whose equation is 

xy = a. 
On this hyperbola connect the points In, -j and (-, nV Let C x be the 




middle point of this chord. It is easy to see that the coordinates of Cx are 
( n u n i) • Then repeat the process with % in place of n, and so on indefinitely. 
An examination of the figure shows at once that the points Ci, C 2 , C 3 • • • 
are moving down the diagonal towards the vertex, but that they never go 
beyond the vertex. Hence the limit of n t is \Ja. This figure further shows that 
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the value of n assumed at the start has no influence on the final result, but 
that if a value of n near ^a be chosen the process will converge much more 
rapidly than if we begin with a value remote from \/a. 

It is interesting to apply this geometrical method to trying to find the 
square root of a negative number. The process does not converge. 

Finally, it should be noticed that this process is exactly that given by 
Newton's method for approximating to the positive root of the equation 



In that method, the curve 



x* - a = 0. 



y = x 2 — a 



is drawn, any value of x is selected, the tangent is drawn at the corresponding 




Flo. 2. 



point of the curve, and its intersection with the .X-axis gives a new value of 
x, with which the process is repeated. In this particular case the equation of 
the tangent at (n, n 2 — a) is 
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y — (n 2 — a) = 2n(x — n) 
or 2nx = y + n? + a 

whence nj = ^(n + - j , if % be the JT-intercept ; 

and the process is repeated. But this is exactly the method described above. 
Newton's method applied to the _p th root gives at once 



(12) 



If we write 



i=)[a>-i)»+s£i]. 



n t 



nj_! — n ( = h it % — V« = €i, -^ = r,- 
it is easy to prove from (12) that 

(P ~ *)3 £ c ; (g ~ l) ff , , 
2n, <€i< 2[n j -( i) -l)^]' 

and r, < ^zl rf_,(l + r i _ 1 )""». 

These formulae are direct generalizations of (8) and (10), for the actual 
and relative error of any approximation . 

Harvard TJkiterbitt, 

Cambridge, Mass. 



